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Abstract. We describe a recent experiment performed with rubidium atoms (^^Rb), aiming at 
studying the coherence properties of a two-dimensional gas of bosonic particles at low temperature. 
We have observed in particular a Berezinskii-Kosterlitz-Thouless (BKT) type crossover in the 
system, using a matter wave heterodyning technique. At low temperatures, the gas is quasi-coherent 
on the length scale set by the system size. As the temperature is increased, the loss of long-range 
coherence coincides with the onset of the proliferation of free vortices, in agreement with the 
microscopic BKT theory. 

1. INTRODUCTION 

The type of phase transitions that a physical system can undergo depends crucially 
on its dimensionality. In a three-dimensional (3D) system, one usually observes the 
emergence of a long range order below a critical temperature. The corresponding order 
parameter can be constant in space, like the macroscopic wave function of a 3D Bose- 
Einstein condensate (BEC) or the magnetization in a 3D ferromagnetic material. In a 
two-dimensional fluid, long-range order is destroyed by thermal fluctuations at any finite 
temperature T . In particular, for a 2D uniform ideal gas of identical bosons it is well 
known that BEC cannot occur at a non-zero temperature. 

However in an interacting 2D Bose gas, there exists a certain critical temperature Tc at 
which a phase transition takes place. Below Tc correlations of the order parameter decay 
only algebraically in space and the system is superfluid. Above Tc the correlations de- 
cay exponentially and superfluidity is lost. The Berezinskii-Kosterlitz-Thouless (BKT) 
paradigm 0,0] associates this transition with the emergence of a topological order em- 
bodied in the pairing of vortices with opposite circulations. For T <Tc vortices are only 
found as part of bound pairs, whereas for T > Tc unbound vortices proliferate and cause 
the exponential decay of the correlations of the order parameter. 

This behavior is characteristic of a variety of 2D phenomena including the superflu- 
idity in liquid helium films the superconducting transition in arrays of Josephson 
junctions [4], and the collision physics of 2D atomic hydrogen [5]. Here we report ex- 
perimental evidence for BKT- type phenomenon in harmonically trapped dilute atomic 
gases. Specifically, we analyze series of interference patterns between independent quasi 
2D gases of ^^Rb atoms. We identify a temperature region close to the predicted T^ 
around which the phase coherence of the gas evolves rapidly with T. For temperatures 
above this transition region, phase defects indeed proliferate. Our observations thus pave 
the way to a direct, quantitative confirmation of the microscopic basis of the BKT theory. 

The paper is organized as follows. We first present (section |2|) our experimental 



configuration in order to set the scale of parameters of our quasi-2D gas. We then 
review the relevant known results for homogeneous (section [Sj and trapped (section |4l) 
Bose gases in two dimensions. In section |5l we present our main experimental results 
concerning the observation of a BKT type crossover in a 2D gas of rubidium atoms. 



2. EXPERIMENTAL SYSTEM 



Our experimental setup has been described in detail in ||^ |7|] (see for other recent 
experiments performed with quasi-2D atomic gases). We start our experiment with a 
magnetically trapped, quantum degenerate 3D cloud of ^^Rb atoms. A ID optical lattice 
with a period of = 3 /im along the vertical direction z is used to split the 3D gas into 
two independent clouds and to compress them into the 2D regime. 

The two clouds form parallel, elongated 2D strips, characterized by the harmonic 
trapping frequencies in the xy plane (Ox/ {2n) = 1 1 Hz and (Oy/{2n) = 130 Hz. The oscil- 
lation frequency in the lattice potential along the z direction is (0^/ {2n) = 3.6 kHz, and 
the tunneling between the two planes is negligible on the time-scale of the experiment. 

The maximal number of atoms in each planar quasi-condensate is ^ 5 x 10^. The 
Thomas-Fermi approximation presented in section|4lyields pc = 10^^ cm^ for the central 
density of each quasi-condensate, and 27?;^ =120 jim and 2Ry = 10 jim for the x and y 
lengths of the strips, respectively. The corresponding chemical potential is ji/h— 1.8 
kHz. Since /i < hco^, the motion along z is frozen out. 

In quasi-2D, the chemical potential can be written as 

where M is the atomic mass, and the dimensionless parameter g can be expressed in 
terms of the 3D scattering length a = 5.2 nm and the extension of the ground state along 
the lattice direction = \^h/ {M(0^ =0.17 /im 



1^78^—^0.15. (2) 



Since g < 1, each planar cloud is a weakly interacting system. The healing length 
= /i/ ^/2^^jI = 1/ \J2gPc — 0.2 /im gives the characteristic size of the core of a vortex, 
and the quantity % j 2g^\^ can be interpreted as the number of atoms in the area 
that are missing in the vortex core. Note that strictly speaking g also depends on the 
relevant energy set by the chemical potential /i |^ 0, O], but for our experimental 
parameters this dependence may be neglected. 



3. THE HOMOGENEOUS 2D GAS 



In this section we give a brief summary of the known results concerning a 2D gas of 
bosonic particles. Consider first an ideal gas of N bosons at temperature T, confined in 



a square box of size L^. Using the Bose-Einstein distribution and assuming a smooth 
variation of the population of the various energy states, we find a relation between the 
spatial density p = N /L^, the thermal wavelength X—h/ {IkMUbTY/^ and the chemical 
potential /i: 

This relation allows one to derive the value of }! for any value of the degeneracy 
parameter pA^. It indicates that no condensation takes place in 2D, contrary to the 3D 
case. In the latter case the relation between p^^^^)[? and }i ceases to admit a solution 
above a critical value of p^^^^A^, which is the signature for BEC. The absence of BEC in 
two dimensions remains valid for an interacting gas with repulsive interactions Ll2. .l3i1. 

In spite of the absence of BEC, the physics of an interacting, homogeneous 2D 
gas is still very rich. One can show that at sufficiently low temperature the gas has 
a superfluid component of density Ps When temperature increases, both the 

thermal wavelength A and the superfluid fraction decrease until the product p^A^ reaches 
the critical value 

r = r, : P.A2=4. (3) 

The superfluid density then suddenly drops to zero, and the gas becomes normal [14]. 
The phase transition also manifests itself in the 1-body correlation function for the 
atomic field g\{r) — (y/'^ir) V^(0)). Below Tc one expects a power law decay: 

T<Tc : gi(r)ocr-V(M') , (4) 

reaching gi (r) oc r~^/^ for T = T^, whereas gi{r) decays exponentially for T > Tc. 

We now discuss the physical phenomenon at the origin of this phase transition [HI, 
0]. We use a simple modeling of the system, which consists in neglecting density 
fluctuations of the gas to deal only with phase fluctuations. At low temperature the phase 
varies smoothly with position and can be expanded in normal modes. This amounts to 
consider only phonon excitations of the system. Using the equipartition theorem, one 
can derive the algebraic decay of the gi function given in (|4|). When the temperature 
increases, other excitations like quantized vortices start to play a significant role. A 
quantized vortex is a point where the density is zero and around which the phase rotates 
by a multiple of 271. For our purpose it is sufficient to consider only singly charged 
vortices, where the phase varies by ±271 around the vortex core. The Berezinskii- 
Kosterlitz-Thouless phase transition, occurring at T = Tc (see ©), corresponds to the 
following mechanism: 

• For T <Tc the free vortices are absent. Vortices exist in the system only in the form 
of bound pairs, formed by two vortices with opposite circulations. The contribution 
of these vortex pairs to the decay of the correlation function gi is negligible, and 
the algebraic decay of gi is dominated by the phonons. 

• For T > Tc, the free vortices form a disordered gas of phase defects, which are 
responsible for the exponential decay of gi. 

Neglecting density fluctuations allows for a qualitative understanding of the BKT 
transition, but does not provide the relation between the total density p and the superfluid 



density p^. This question has been addressed in fi?] in the hmit of an ultra weakly 
interacting gas (ln(ln(l/g)) 1). The case of weak, but more realistic interactions, 
has been addressed numerically in [16] and gives Ps/p = 4/ln(^/|) 0.5 for our 
parameters (the dimensionless number ^ is found numerically 380). 



4. THE TRAPPED GAS 

Consider now a gas of N particles confined by the potential V (r) = M( a)^x^ + (Oyy^) /2 in 
the xy plane. The presence of a trap modifies the density of states and BEC is predicted 
to occur for an ideal gas when the temperature is below the critical temperature Tq [17]: 

However it should be pointed out that condensation is a very fragile phenomenon in a 
2D harmonic potential. To illustrate this we calculate the spatial density p (r) using the 
semiclassical approximation. Starting with the expression for the phase space density 

.... 1 1 

w{r,p) = 



^'exp((^+y(r)-M)/W-l' 
we obtain 

p(r) = J wir.p) d^p = -i^ln(l . 

Taking ^ to reach the condensation threshold, we find 

Pmax(r)A2 = -ln(l-^'-^(^)/^«^o). 

Integrating Pmax(^) over the whole space, we recover the threshold ©. However, we 
also note that Pmax(O) = ^. This suggests that condensation in a 2D harmonic poten- 
tial requires an infinite 2D spatial density at the trap center, a criterion which cannot 
be fulfilled as soon as repulsive interactions play a significant role. Note that this is in 
sharp contrast with the case of a 3D harmonically trapped Bose gas, where condensation 

occurs when Pm^^ (0):^2.612/A^, which is finite and (within the semi-classical approx- 
imation) equal to the threshold density in a homogenous system. We refer the reader to 
for a detailed discussion of this subject and a collection of recent references. 
In the limit of large atom numbers and low temperature, we can still determine the 
equilibrium shape of the gas by looking at the balance between the trapping potential 
and the repulsive interatomic potential (Thomas-Fermi approximation). We assume that 
the atomic gas forms a quasi-condensate aizHIIESl, so that ((f t(r))2(v/(r))2) ~ 

((i//"'*(r)i/lr(r))) = p^(^). We then get an atomic distribution varying as an inverted 
parabola 

Pc(r) = Pc(0) I 1 - ^ - 4 I ^8 PciO) = M 




where the chemical potential /i and the radii of the clouds 7?/, i = x,y are: 



In the regime where Ng ^ 1, the chemical potential is much larger than h(b and 
the dynamics of the gas at low energy is dominated by phonons. One then recovers 
several features of the homogeneous 2D gas [23]. In particular for <^ r ^ Rx.Ry, the 
correlation function gi{r) = {y/^{r) V^(0)) decays algebraically at low temperature with 
exponent l/(pc(0)A^) (see eq. (2.59) in [23]). A crossover to a true BEC occurs when 
the decay of gi over the sizes Rx^Ry of the quasi-condensate is negligible. 

Harmonically trapped gaseous samples are finite in size and have non-uniform den- 
sity. For these reasons, one does not expect to observe a sharp BKT transition to superflu- 
idity, but rather a smoothed out crossover. However, an interesting question is whether, 
at least in principle, in the thermodynamic limit one could observe a signature of a 
sharp transition in a harmonic trap. For simplicity we consider an isotropic trapping 
(cOx = cOy = cb) and define the thermodynamic limit by taking N ^ oo and 0) ^ 0, while 
keeping Nco^ constant. We can then duplicate the Kosterlitz-Thouless argument |2] to 
find the probability for having a free vortex appearing close to the center of the quasi- 
condensate under the influence of thermal fluctuations. We choose as a sample region 
the central disk of radius eR, with £ <C 1 and kept constant as 7? ^ oo. The kinetic energy 
of the gas due to a free vortex in the center of the trap is related to the superfluid density 
Ps{r), which we assume to be close to the quasi-condensate density Pc(^) [Z^ 

The entropy associated with the possible choices for the location of a free vortex of size 
^ in the disk of radius £RisS = Ub ln((£7?)^/(^^) = IUb [ln(7?/(^) +ln(£)]. Therefore in 
the limit R/t,^\/£^\,i\it free energy F — Ek — TS associated with a free vortex in 
the central disk of radius eR varies as 

^ = \{psm^-A)HR/^). (6) 

We thus expect a rapid variation of the average number of vortices present in the central 
region when p^(0)A^ varies around 4. When F > (i.e. p^(0)A^ > 4) the probability to 
get a free vortex in this region is small, as the energy cost for a free vortex dominates over 
the gain in entropy. On the contrary for F < 0, i.e. p^(0)A^ < 4, we expect a proliferation 
of free vortices around the center of the trap and a jump of the local superfluid density to 
zero, in analogy with the homogenous system. Finally, note that even for p,(0)A^ > 4, 
free vortices can appear in the periphery of the condensate, because the energy cost for a 
off-centered vortex is lower. The role of these outer vortices on the coherence properties 
at the center of the condensate remains to be studied. 



FIGURE 1. Examples of interference patterns between two quasi 2D Bose gases. Left: low temperature 
pattern; Middle: high temperature pattern; Right: example of a dislocated pattern revealing the presence 
of a vortex in one of the two planes. 

5. EXPERIMENTAL STUDY 
5.1. Matter wave heterodyning 

In order to study the phase coherence of our quasi-2D gases, we have implemented a 
matter- wave heterodyning technique. When all confining potentials are suddenly turned 
off the two clouds expand predominantly perpendicular to the xy plane. As they overlap, 
a 3D matter wave interference pattern forms with a spatial period D — Inht / {md) along 
the z direction, where t is the time-of-flight (TOF) duration ll25ll . After t = 20 ms, the 
projection of the 3D interference pattern onto the xz plane is recorded on a camera, using 
a resonant probe laser directed along y. 

Three examples of interference patterns are shown in fig.E The left pattern is typical 
of what is obtained at very low temperature, with nearly straight fringes. The middle 
pattern has been obtained for a larger temperature and the presence of phase fluctuations 
in the two planes is revealed by the waviness of the interference fringes. Finally we 
occasionally observe patterns with one (or several) dislocation(s), such as the one shown 
in the right of fig. [H We interpret such a dislocation as the signature of a free vortex in 
one of the two atomic clouds The matter wave heterodyning technique thus 

offers a way to study both long wavelength excitations (phonons) by analyzing the 
smooth variations of the fringes of the interference patterns, and the point-like defects 
such as free vortices by looking for sharp dislocations in the fringes. 



5.2. Extracting gi from interference patterns 

We use in this paragraph a simple modeling, where we assume that the two gases 
have the same uniform amplitude i/Zq and independently fluctuating phases (pa{x^y) and 
%{x,y). The interference signal recorded after TOF is 



(7) 



with ^ 

Ly J-Ly/2 

At this stage the integration distance Ly is arbitrary, and it can differ from the total length 
2Ry along the y direction. For a quantitative analysis of long-range correlations as a 
function of temperature, we adopt the method proposed in ll28ll . The idea is to integrate 
the coefficient c{x) appearing in ^ over a variable length L^: 

1 rLx/2 
C{L^) = — / c{x) dx 

Lx J-Ljc/2 

and average \C{Lx) ^ over many images recorded in the same conditions. Using the fact 
that the phases % and % are uncorrelated, we obtain: 

(|C(L,)|2) = ^Jl{c{x)c*{x'))dxdx' 

(^gK%{x,y)-%{x',y')^ ^gi{n{^y)-cpb{x,y)^ dxdx'dydy' 

fLx/2 / I 

■ / \gl{x,0)\^dx oc - 



1 /-ws , /n2« 

, (9) 



L 



where we have assumed Lx^ Ly. The long-range physics is then captured in a single 
parameter, the exponent a. It is straightforward to understand the expected values of a 
in some simple cases. In a system with true long-range order, g\ would be constant and 
the interference fringes would be perfectly straight. In this case a = 0, corresponding 
to no decay of the contrast upon integration. In the low temperature regime, where g\ 
decays algebraically (see ©) the exponent a coincides with that of g\. In the high 
temperature case, where gi decays exponentially on a length scale much shorter than Lx, 
the integral in © is independent of Lx. In this case a = 0.5, corresponding to adding 
up local interference fringes with random phases. The BKT mechanism corresponds to 
a transition between a power law with exponent (l/{ps?i^) < 0.25) to an exponential 
decay of gi. It should thus manifest itself in a uniform system as a sudden jump of a 
from 0.25 to 0.5 when the temperature varies around Tc. 



5.3. Experimental results 

In order to apply the method of [28] to our data, we fit the density distribution recorded 
after TOP with a function 

F{x,z) = G{x,z)[l+c{x) cos(27rz/D + 0(x))] , 

where G(x,z) is a gaussian envelope, and we proceed with the experimental c{x) as 
discussed above. We average \C{Lx)\^ over ^ 100 images and fit it with the function 
L~^^. We have plotted in fig. |2K the results for a, as a function of the average central 
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FIGURE 2. (a) Fitting exponent a as a function of the average central contrast cq. (b) Fraction of 
images showing at least one dislocation as a function of cq. 



contrast cq of the images. We choose the latter quantity as a temperature label since 
we feel that it is less model-dependent than other estimates, based for example on 
the tails of the atomic TOF distribution. Starting at high temperatures, for values of 
Co ranging from a few % (experimental limit of detectability) up to about 13%, a is 
approximately constant and close to 0.5. When the temperature is reduced further, a 
rapidly drops to about 0.25, and for even lower temperatures (larger cq) it levels off. We 
thus clearly observe a transition between two qualitatively different regimes at high and 
low temperatures. The values of a above and below the transition are in agreement with 
the theoretically expected jump at the BKT transition in a uniform system. 

Figure |2b shows the number of sharp dislocations at different temperatures. For 
the count we consider only the central, 30 jim wide region of each image, which is 
smaller than the length of our smallest quasi-condensates. We observe a clear onset of 
dislocation proliferation with increasing temperature, and this onset coincides with the 
loss of quasi-long range coherence described by the rapid variation of a. 

We can also independently estimate the cloud's temperature and quasi-condensate 
density at the onset of quasi-long-range coherence. For images with cq = 0.15, we get 
Pr?i^ = 6 ± 2 [7]. This is in agreement with the BKT theory, and we refer the reader to 
infll for a further discussion on the difference between pc and p^. 



5.4. Discussion and summary 

The experimental procedure presented here constitutes a powerful technique to study 
the properties of quasi-2D atomic gases at ultra-low temperature. Our results show that, 
in agreement with the prediction of the BKT theory, there exists a relatively narrow 
temperature domain over which the coherence of our trapped planar Bose gas varies 



rapidly, with the exponent a in (|9l) switching from ?^ 0.25 (low temperature side) to ^ 
0.5 (high temperature side). In addition the matter wave heterodyning technique allows 
for the direct observation of dislocations of the interference pattern due to vortices, and 
it supports the notion that the unbinding of vortex-antivortex pairs is the microscopic 
mechanism destroying the quasi-long range coherence in 2D systems. 

However, we should point out that the quantitative agreement between our results 
and the predictions of [28] derived from BKT theory might be partly fortuitous. Indeed 
though we concentrated on the quasi-uniform part of the images, the geometry effects in 
our elongated samples (R^ ^ l2Ry) could still be important. In particular the integration 
length Ly in our experiment is equal to the full width of the system 2Ry in that direction. 
The condition ^ Ly used to derive ^ then implies that we probe the planar system 
on a distance scale that is large compared to its other dimension 2Ry. Over this length 
scale the system acquires some quasi-one dimensional features and gi is expected to 
decay exponentially, possibly with a very large coherence length (see e.g. [29]). If this 
is indeed the case, our experiment is detecting a rapid change of the decay length of the 
exponential, rather than a change in the functional form of g\ given in (|4|). The elongated 
nature of our samples also changes the quantitative calculation of the threshold for the 
apparition of a vortex (see section |4|). Finally we note that we detect here only a subset 
of vortices, i.e. those which are well isolated and close to the center of the cloud. The 
heterodyning technique is not well adapted to the detection of large numbers of vortices, 
that are expected for temperatures notably above Tc. In the latter case indeed, the fringe 
pattern becomes hardly detectable and thermally activated phonon modes with a very 
short wavelength along x can also have a significant contribution. Our work should then 
be considered as a first step in the investigation of the BKT mechanism in planar atomic 
gases. A quantitative theoretical analysis of our system will probably require a numerical 
study, for example along the lines of [30]. In parallel, future experimental studies could 
be performed with more isotropic samples, and could use a "slicing method" similar to 
ll25[l in order to extract local information from the center of the cloud. 
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